SCALING LIMITS OF RANDOM SKEW PLANE PARTITIONS WITH 
ARBITRARILY SLOPED BACK WALLS 



SEVAK MKRTCHYAN 

Abstract. The paper studies scaling limits of random skew plane partitions confined 
to a box when the inner shapes converge uniformly to a piecewise linear function V of 
arbitrary slopes in [—1, 1]. It is shown that the correlation kernels in the bulk are given 
by the incomplete Beta kernel, as expected. As a consequence it is established that the 
local correlation functions in the scaling limit do not depend on the particular sequence of 
discrete inner shapes that converge to V. A detailed analysis of the correlation kernels at 
the top of the limit shape, and of the frozen boundary is given. It is shown that depending 
on the slope of the linear section of the back wall, the system exhibits behavior observed 
in either |OR07| or |BMRT10| . 
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1. Introduction 

1.1. Background. 

1.1.1. Skew plane partitions. Given a partition A, a skew plane partition with boundary 
A confined to a c x box is an array of nonnegative integers vr = {tTj ,,■} defined for all 
l<i<c, l<j<d, {i,j) 4- -^j which are non-increasing in % and j. One way to visualize 
this is to draw a c x (i rectangular grid, remove the partition A from a corner, and stack 
TTjj identical cubes at position (i, j), as in Figure [T] The number of cubes |7r| := X^^jj is 
the volume of the skew plane partition tt. 



Figure 1. A skew plane partition and the corresponding partition A. Here 
A = {4, 4, 3, 3, 3,1}. 




From Figure [T] it is easy to see that skew plane partitions can be identified with tilings 
of a certain region of with 3 types of rhombi (see |OR03j or |OR07j for details and 
other correspondences). Scale the axes in such a way that the centers of horizontal tiles 
are on the lattice Z x The letter t will be used for the horizontal and h for the vertical 
coordinate axes in this plane. For a partition A, let h\{t) encode the boundary of A. More 
precisely, let Mi < M2 < . . . < Un-\ denote the t coordinates of the corners on the outer 
boundary of the Young diagram A as shown in Figure ll Define hx{t) as 




n-l 

6a (t) := t + 2 Y,^-\)\t - u,)9{t - M,), 

i=l 
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where 9(t) is the step function 

e{t) 



1, t>o 

0, t < 



We will call bx{t) the back wall. Notice that 6A(^is a piecewise linear function with slopes 
in [—1, 1]. The graph of h = ^bx{t) (see Figure^ gives the inner shape of the skew plane 
partition from Figure [T] 



Figure 3. The graph of ^bx{t) when A = {4, 4, 3, 3, 3, 1}. 

h 




1.1.2. The thermodynamical limit; cases studied before. For q G (0, 1) introduce a prob- 
ability measure on skew plane partitions with boundary A and confined to a c x c/ box 
by 

Prob{7r) oc g'^L 

Given a subset U = {{ti, hi), . . . , (t„, hn)} C Z x |Z, define the corresponding local 
correlation function px q{U) as the probability for a random tiling taken from the above 
probability space to have horizontal tiles centered at all positions (tj, K)^^^. 

Okounkov and Reshetikhin studied the thermodynamical limit q = e"'" — )■ 1 of this 
system in several cases. In |OR03j they studied the case when A is the empty partition 
and the size of the box is infinite (i.e. plane partitions are not confined to a finite c x d 
box). 

In |UR07j they showed that in general, for arbitrary A, the correlation functions are 
determinants. 

Theorem 1.1 (Theorem 2, part 3 |OR07| ). The correlation functions p\^q are determi- 
nants 

P\,q{U) = det{Kx,q{{ti, hi), {tj, /ij)))l<ij<n, 

where the correlation kernel Kx^g is given by the double integral 

(1) Kx,g{{tl,hi),{t2,h2)) = 

1 [ [ '^6,(2:, ti) ^_^,^^,^(,^^_^^^^^^_l_,^(t^)^dzdw ^ 

(2vri)2 J,^c, Jw&c^ {w, t2)z-w zw ' 
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where h\{t) is the function giving the hack wall corresponding to A as in Figure^ 

^-,h^{z,t) 



'^+,h,{z,ty 



2 

— 1 „— m\ 



= n (1 

m<t,m(iD~ ,mSZ+| 

m G means the hack wall at t = m, i.e. bx{t) at t = m has slope =f1, and Cz (resp. 
Cw) is a simple positively oriented contour around such that its interior contains none 
of the poles o/$fe^(-,ti) (resp. of ^b^{-,t2)~^ ). Moreover, if ti <t2, then Cz is contained 
in the interior of C^, and otherwise, C^ is contained in the interior of Cz- 



Using Theorem 1.1 they studied the hmit g — )■ 1 of the statistical system when the 
number of corners in the partition A stays finite and showed that it exhibits the hmit 
shape phenomenon. More precisely, let {r^j^j^ be a positive sequence which converges 
to 0. Set qk = e"*"*. Let Xk(t) be a collection of partitions confined to Ck x dk boxes with 
corners at positions < < • • • < u^ti, and let bxi^{t) be the corresponding functions 
giving the back walls. Note, that the positions of corners depend on k but the number 
of corners does not. Consider the probability measure on skew plane partitions vr'^ with 
boundary Afc and confined to a Cfc x dk box by 

(3) Pro6(7r^) ocgf''. 

The characteristic scale is r^^, so the system should be scaled by r^. Let Bx^t) 
denote the function corresponding to hx^t) when scaled by r^, i.e. Bx^t) := rkbx^.{^)- 
Given constants Vq < Vi < . . . < Vn-i < Vn such that Vq < < Vn, Okounkov and 
Reshetikhin studied the limit k ^ oo when — )■ 0, VkU^'' — )• Vi for all 1 < i < n — 1, 
-^fkCk — Vq, and -^Tkdk — )■ Vn- In other words, in the scaling limit the back wall is a 
piecewise linear curve with line segments Vi < r < V^+i of slopes ±1 (see Figure lij). 




They showed that the scaled height function r^.Tr'^ of a random skew plane partition 
converges to a deterministic shape. In the two dimensional formulation of the system 
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as random tilings of the plane by lozenges two different types of regions appear: those 
where only one type of tiles appear, called frozen regions, and those, where tiles of each 
type appear with positive density, called liquid regions. They showed that the correlation 
functions in the limit are described by a determinantal process given by the incomplete 
Beta kernel in the liquid region. They gave a characterization of the boundary between 
the frozen and liquid phases that arise and showed that in the limit the correlation kernel 
K is given by the Airy kernel at generic points of the frozen boundary and by the Piercey 
process at the singular points of the frozen boundary, which appear near the outer corners 
(see Figure |5]). 

Figure 5. The frozen boundary in |OR07j . 




While in |OR03j and |OR07] the limiting back wall consisted of a piecewise linear 
function of slopes ±1, in |BMRT10] we studied the case when the limiting back wall is 
piecewise periodic in such a way, that in the scaling limit it converges to a continuous 
piecewise linear function of rational slopes strictly between —1 and 1 (see Figure |6]). We 
showed that the correlation functions in the limit are given by the same processes as 
before, both in the bulk and at generic points of the frozen boundary. Unlike the |UR07j 

Figure 6. The periodic back wall in |BMRTinj . 




case, the frozen boundary doesn't develop cusps (see Figure [7]); hence, the Piercey process 
does not appear. Another difference is that in this setup the liquid region extends to oo 
everywhere on the back wall (see Figure [T]), and the local statistics is given by the bead 
process of Boutillier |Bou09j when you move high up on the wall. The bead process also 
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appears high up, near the corners, with various parameters depending on the slope of the 
approach to the corner. 



Figure 7. The frozen boundary in [BMRTlOj . 




1.2. Main results. The present work studies the general case when the limiting back 
wall is an arbitrary continuous piecewise linear function of slopes in [—1,1]. This is a 
generalization of |OR07] where slopes were ±1 (lattice slopes) and of |BMRT10] where 
slopes were in (—1, 1) (non-lattice slopes). 



Let rk,qk,Ck,dk, Xk,bxf^ and Bx^ be defined as in Section 1.1.2, In Section 1.1.2 
the number of corners in was independent of k, however here there is no restric- 
tion. Consider the limit — )■ oo when — )■ 0, r^Ck — )■ C, Vkdk — )■ D, and the curves 
B\^{t) converge point-wise and uniformly to a continuous piecewise linear function V{t) 
with slopes in [—1,1]. The following notation will be used throughout the paper. Let 
Vo < . . . < Vj^-i < T <Vj^ < . . . <Vn and . . . , /3n G [—1, 1] be such that V{t) is linear 
for r G \Vi-i, Vi\ with slope Note that jr will denote the index where Vj^-i < t <Vj^. 
In order to simplify notation, from now on the index r from jV will be dropped. For 
convenience also define (3q = —1 and f3n+i = 1, and require that /3j 7^ for all i. 



Figure 8. A frozen boundary. 




In such a scaling limit the statistical system of random skew plane partitions dis- 
tributed according to (Isj) behaves as follows: 
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(1) The system exhibits the hmit shape phenomenon. 

(2) Near sections of the back wall where the slope is ±1 the disordered region is 
bounded above as in |OR07] . Near sections with slope in (—1, 1) it grows infinitely 
high as in [BMRTinj (see Figure g). 

(3) The number of connected components of the frozen boundary is one more than 
the number of outer corners where at least one of the slopes at the corner is a 
lattice slope. Note: outer corners are those where V'{t—) < V'{t+). 

(4) The frozen boundary develops a cusp for each such outer corner. 

(5) The correlation functions are given by determinants with the incomplete beta 
kernel in the bulk, the Airy kernel on the frozen boundary, and the Piercey kernel 
near the cusps, as in |OR07] . 

(6) High up, near the sections of the back wall which have non-lattice slopes, the 
correlation functions converge to the bead process as in |BMRT10] . High up, near 
the sections of the back wall which have lattice slopes, the region is frozen. 

(7) When approaching a corner and simultaneously going up along the frozen bound- 
ary, depending on the relative speed of the approach and the slopes of the back 
wall at the corner, either a frozen region or the bead process can be seen, with the 
density of the beads depending on the relative speed of the approach. 

In this case, unlike |UR03j . jnR06] . |UR07] and |KU07] . the boundary of the limit 
shape is not an algebraic curve. 

Cusps on the frozen boundary can exhibit previously unknown behavior, which will 
be studied in future articles. One such example is when cusps approach each other as 
shown in Figure |9| 

Figure 9. The frozen boundary when the back wall is given by = 
-12.1, -12, -8.1, -8, -4, and /3 = -0.9, -1, -0.9, 1, -1. The left picture 
is the right picture zoomed in inside the circle. 




It is also shown in this paper that the local statistics in the scaling limit is independent 
of the intermediate approximations. Let and be as before. Let t'l G Z, t2 ^ ^) 
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G ^Z, and G be sequences such that 
(4) hm t'l = hm t2 = T and hm h'l = hm = X- 

k—^oo k—^oo k—^oo k—^oo 

The following theorem is proven. 

Theorem 1.2. Let Xk{t) and fik{t) be two sequences of partitions with corresponding 
scaled boundary functions Bx^^ (r) and (r) . // in the limit k ^ oo both B\^, (r) and 
A\^{t) converge uniformly to the same piecewise linear function V{t) corresponding to 
the back wall, then 

fe— ^-oo fe— ^-oo 

In other words, limfc_5.oo Kx^,q^{{t\, h\), (tg, ^2)) only depends on the limiting back wall 
V{t). It does not depend on the intermediate steps Afc(t). In particular this is true for 
the one point function hmfc_,oo PAfe,<?fc(^i, hf) = limfc^oo Kx,,,q,^{{t\, h^), {tf, h^)), which gives 
the slope of the limit shape at the point (r, x). Thus, it is obtained that the limit shape 
is also independent of the intermediate steps. 

1.3. Outline of the structure of the paper. One of the main steps in what follows 
is to understand the asymptotics of In $6^^ . This is done in Appendix lAl The reason for 
separating this from the rest of the paper is that the argument actually holds in greater 
generality; one can consider back walls where V{t) is an arbitrary continuous Lipshitz 
function with constant 1. 

Section [2] studies the critical points of the asymptotically leading term in the integral 
formula giving the correlation kernel. The main result is that the number of non-real 
complex critical points is or 2. The general strategy for proving this is the same as that 
from |BMRT10] : obtain the number of non-real complex critical points when x — )■ 00 and 
from this deduce the result for finite x- However, the argument in |BMRT10] could not 
be used here, since it heavily relies on the frozen boundary having a certain simple shape 
with only one connected component (see Figure [?]). 

The results of Section |2] are used in Section |3] to apply the saddle point method and 
obtain the limit of the correlation kernel in the bulk. 

Section |4] gives a detailed study of the frozen boundary. Such a study has not been 
carried out before. 

1.4. Acknowledgements. I am very grateful to Peter Tingley for a suggestion on which 
the proof in Appendix |X] is based. I am very grateful to Cedric Boutillier and to an 
anonymous referee for many suggestions to improve the presentation of this paper. I 
am very grateful to Nicolai Reshetikhin for his guidance. I am also very grateful to 
Alexei Borodin, Cedric Boutilier and Peter Tingley for many useful discussions on the 
subject. Lastly, I would like to thank the organizers of the Park City Mathematics 
Institute Summer School on statistical mechanics, in July 2007, where much of Appendix 
[Al was written. 
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2. Critical points of the asymptotically leading term in the integral 
formula giving the correlation kernel 



2.1. The function S{z). Lemma A.l gives an asymptotical formula for {z, t) defined 
in (|2]), when the partitions are not restricted to a box. If they are bounded to a box in 
such a way that r G [Vq, Vn], then from the definition of $6^^ {z, t) it is easy to see that if 

bxf, {t) are defined by the conditions 

bxi,(t) is continuous, 

K(t)= -1, nt<Vo, 
Kit)= Kit), Vo<nt<Vn, 



then $^ {z,t) = {z,t). Moreover, in the limit — )■ 0, the sequence Bx^{t) (i.e. the 
scaled bxi^{t)) converges to the function V{t) defined by the conditions 

V{t) is continuous, 

V'{r) = -1, T< Vo, 

V(r)= V{t), Vo<r<Vn, 

V'{r) = 1, Vn < T. 

Let t^, h^, i G {1, 2} be as in (|4]). The previous two statements, together with Corollary 



A. 2 give that for plane partitions restricted to a box in such a way that r G [Vq, Vn] 



lim rfcln$fe = lim rfcln<l>r {z,t^) = 

k—^oo " k—^oo k 



- [ --(l + \/'(M))ln(l-e*Vi)dM+ /" " --(1 - F'(M)) In (l - e'^^z) dM. 
Jvo 2 2 

Now, from the integral formula ([!]) conclude that in the limit — )• 0, r^t^' — )■ r, 
fkhf — the correlation kernel can be written as 

(5) Kx,,,M,h1),{4,hl)) = —— / / e ^ 



(27ri)2 J J z — w 

where the contours of integration are as in ([T]), and Sr^xi^) is defined as 

(6) Srjz):=- r -1(1 + V'iM))\n{l- e^'z-')dM 
JVo ^ 

+ ^ " -^(1 - y\M)) In (1 - e-'^'z) dM - ln(z) [x - ^l^(r)^ . 



2.2. Number of complex critical points of S{z). In Section|3]the asymptotics of the 
correlation kernel will be studied by the saddle point method. Since the asymptotically 
leading term of the correlation kernel is given by the function Sr^x^z), to use the saddle 
point method the critical points of Sr^^^z) need to be studied. The goal of this section is 
to show that S^-^xi^) exactly 2 or non-real complex critical points. The number of 
such critical points depends on the position of the point (r, x). 
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To simplify formulas, from now on the subscripts in Sr^xi^) will be omitted unless 
that might cause ambiguities. 

2.2.1. Formulas for zS' and z{zS'y . In the study of critical points of S{z) formulas for 
z'^^^j^ and z-^{z'^^^^) will be needed. The reason for working with these functions instead 
of S' and 5"' is that expressions for S' and S" are more complicated. 

It follows from ^ that 

- (X- - \v{t)) 

- 1 + ft) (I'r^^) - i(i + ft) i.> (r' - 



+2(^-^^)^HT3^j + 5: 2^1-/5^)1^1^^) 

\ / i=j+i ^ / 

-(x-^^r)). 

Using Viyi) = Viyi^i) + Piiyi — Vi_i), ^^^^ can be rewritten as follows: 
^ci5(z) _ , / -ze-^^ - 1 ^^ 1/1 , « / -26"" - 1 



1 / 7p-^^ - 1 \ 1 / 
V-(l + /3i)ln --(l + /3,)ln 



From here it is easy to obtain 



4z (^l^w) = ^(i+ft)F^+i:^(ft.«-ft)F^4'^-^»'i^-7^- 

2.2.2. Critical points of S{z) are away from [e^'-^, e^'] z//?, 7^ ±1. If -Zcr is a critical point 
of ^(z), then 2cr^^^ = 0. 

Lemma 2.1. //i 7^ j zs such that /3j 7^ ±1, then there does not exist zq G [e^*-i,e^'] snc/i 
i/iai lim^^^o^;^^ = 0. 

Proof. Suppose Zq = x + ie, where x G (e^'-^e^*) and \e\ ^ 1. For such zq all 
terms in ([T]), except perhaps one, have imaginary parts close to zero. More precisely, 

^ / dSjzo) ( z^e-^^- \ W 

^r^-^l" Uoe-^--l JJ^'^^'^' 
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where c is the coefficient of ln( ^^°f v-^^ in Since In is the branch of the logarithm 
with an imaginary part in (— vr, vr) with a cut along ]R_, then 

^oe-^^ - 1 



53 cln 



ZqC-^^-^ - 1 



±c7r + 0(e) 



which in turn implies ^^(^o^^^jf^) = ±c7r + 0{e). From ([T]) it follows that c = ^(1 + j3i) 
or c = |(1 — f3i). If f3i 7^ ±1, then c 7^ 0. Thus, there is no zq G (e^'-^e^*) such that 
0. The only remaining points are Zq = e^'-'^ and Zq = e^', but these are 



lim. 



•-z^zo ^ dz 

singular points of z'^^^^ and therefore Ymiz^zo z'^^^j^ ^ 0. 



□ 



Notice, that if /3j ■ 
In this case ^[zq'^^^'^' 



dz 



1 and z > j, or /3j = — 1 and i < j, the coefficient c is zero. 
= 0(£:), and there may be critical points of S{z) in the interval 



2.2.3. Critical points when % — > cxd. 

Lemma 2.2. Fix r G (V^-i, V^)- -^or sufficiently large x, St,x{^) has no non-real complex 
critical points if Pj = ±1 and exactly two non-real complex critical points (which will be 
complex conjugates) if f3j 7^ ±1. 



Let us ffist prove the following lemma, which will be used in the proof of Lemma 2.2 



Lemma 2.3. Let {xi}™^, {yj}S=\, {fi}i='i ^ \ m 7^ 6e real numbers such 

that 

Xm, < ■ ■ ■ < X2 < Xi < w < yi < y2 < ■ ■ ■ < ym2- 
Let C be a complex number with '^C, > 0. Define 

/ii = angle(C - Xj+i, C - Xi), \fi = 1, 2, . 
Ui = angle(C -yuC- Vi+i), = 1,2,. 

and 

a = angle(C — Xi,( — "w). 



,m2 ■ 



L 



(See Figure 10 for an illustration of the setup on the complex plane.) There exists e > 



Figure 10. Setup of Lemma 2.3 

C 

/^mi-l 




... X2 xi uo yi y2 ' ' ' V'^^ 
such that if — w\ < e, then 

mi — 1 r7i2 — 1 

(9) ^ ^iHi + ma + ^ OiZ/i = 



i=l 



i=l 
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if and only z/ G M and ( > w. 

The same holds if a is replaced by a := angle(C — wX — Vi) o^'^ ( > w is replaced by 
C<w. 

Proof. Let e = yi — w. U ( G M., ( > w and |^ — w| < e, then w < ( < yi. Thus, 
all the angles fXi, z/j and a are zero, and ^ is true. 

Let us prove the converse. Suppose $e > such that (|9| implies C £ I^j C > If 
^ G M and Xi < ( < w, then all the angles /Xj and Ui are zero, but a = n and ^ cannot 
hold. It follows, that there must be a sequence Q e C\M such that 

(10) lim Ci=w 

and ([9]) holds for ( = Ci^Vi Notice that the angles a,Hi, and z/j depend on ( and that 

(10) implies 

(11) lim fii{Ci) = and lim uilQ) = for all I. 

These, together with the assumption that ^ holds for C = '^h imply 

(12) lim a(0) = 0. 



Define ri{Q := angle(C — wX~ ^C) and 7(C) := angle(C — C — yi) (see Figure 11 ) 
Using this notation. 



Figure 11. 




X2 Xi 



^ J^C yi y2 



tan(a(Ci)) 





lb 


— Xi 


cos(r/(0)) 


10- 


w 


+ 


w 


- Xi 


sin(?/(Ci)) 



Writing a similar expression for a(Ci) + j^iid), and using (11) and (12), obtain 
/ii(Ci) + a(Ci) ,. tan(/ii(Ci) + a(Ci)) 



(13) lim 



a(0 



lim 

j-s>oo tan(a;(Ci)) 

|w-X2| cos(??(Ct)) 



lim 



|Ci-u'|+|w-X2| siii(}?(Ci)) 

[w-xil cos(??(C»)) 
\(i-w\+\w-xi\sin{rj{Ci)) 



lim 



3 + sin(r/(0)) 



|«)-X'2| 



It follows from 1^ and that 
(14) lim 7/(0) = lim 7(0 



TT 

2' 
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Thus, (13) gives that hm, 

easy to see that u 
/ = 1, 2, 3, ... , mi. 



Mi(C»)+a(C») 



1, and hence that hm. 



Mi(C») 



easy to see that using the same argument it can be shown that hmj^oo 



0. It is 
for all 



It follows from ([TO]) and ([14]) that 

^i(C) 



tan(i/i(Ci)) ^ |^,_j,,|+|j,,_j,2|sin{7(0)) 



ll/l - 


i/2 






+ 


Z/1 


-?/2| 



tan(a(Ci)) 



lim 



|C,-«)| + |ui-a;i I sin(r){0)) 



cos(7(Ci)) _ 




-1/2I 




cos(r/(0)) 


- y\\ 


+ bl 


-1/2I 












\w 


- V\\ 


+ \y\ 



lim 



sm 



-7(0)) 



sin(|-r/(0)) 



lim 



Similarly, lim 



10 - ^^1 

\v\ - V2 \ i^'^ 10 - yi\ 

for all / = 1, 2, 3, ... , 1712- Combining the results gives 



lim 



1=1 



MO 

a(o; 



"12-1 / > \ 



which is a contradiction to the assumption that ^ is satisfied for C = fo^' ^• 

This proves the first statement in the lemma. The second statement follows by sym- 
metry. □ 



Proof. [Proof of Lemma 2.2 Since z 0, studying the critical points of S{z) is 
equivalent to studying the solutions to z-^S{z) = 0. The real part of the equation, i.e. 
= 'Si{z^S{z)) implies that if x is very large, then z must be very close to e'^ or to e^' 
for some /. 

The imaginary part = '^{z-^S{z)) is equivalent to 



1 

(15) = -^-(l + /3Jangle(z-e 



- ^(1 + /3j)angle(z - e^'-\z - e") 



1 1 

+ -(1 - /3j)angle(2 -e\z-e^^)+ ^ -(1 - /3^)angle(2; - e^" 



)• 



m=j+l 



For arbitrary real numbers x, a, 6, if x ^ [a, 6], it is immediate that 

lim angle(2; — e'^, 2; — e^) = 0. 



This implies that if z is near e then all but the two terms in (15) where e ™ appears 
are close to zero. The sum of the angles in the remaining two terms is vr, and both have 
coefficients of the same sign. If none of those two coefficients is zero, then if z is sufficiently 
close to e^"", the RHS of (15) cannot be zero. 



If one of the two coefficients is zero, then in order for (15) to hold, z must be real. This 
follows from Lemma 2.3, For example, in the case Pm 7^ — 1? Pm+i = — 1 and m < j — 1, 
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setting C = z, 
(16) 
and 
(17) 



Xi 



w 



yi 



yj—m—l yj—m yj—m+1 





. . . e^™-! 




f 1 


J-mi — 1 • • • 




to 


t)l 


ll/3i 






l + /9m + 2 


2 


2 


2 


2 



ym2 



lift- 1-k 



• • • t)m2-l 



2 ' 



Lemma 2.3 gives that (15) implies z must be real. 



If z is near and /3j = ±1, then exactly one of the coefficients of angles in (15) 
containing 2; — e"^ is zero, and again it follows from Lemma 2.3, with the parameters set 



up similarly to (16) and (17) but this time with w = e^, that z must be real. 



If z is near e'^ and Pj 7^ ±1, then using the methods from the proof of Proposition 
3.1 of |BMRT10] . it can be shown that there are two possibilities for z, and these two 
complex conjugate critical points of S{z) can be asymptotically calculated. Formulas for 



the critical points are given in Lemma 2.4 
Lemma 2.4. In the limit x — t- 00 



□ 



(a) Ifr G (V^-i, Vj), f3j 7^ ±1 is fixed, then the asymptotics of the non-real complex critical 



points is given by Zc 



', where 



^ = e-Xe±i4(i+/3.)-Q 

j=0 



2sinh 



T-V, 



(h) If r = Vj-i + 5, X ^ 00 and 6 in such a way that 

p ^ ex-x(^-i)|^|i-|(ft-/3i-i) 

is fixed, with 



^(i-i) 



Y[ |2sinh( 



i=0 



then the critical points behave as z^r = where s is a solution to the equation 

(18) „ _ ,±i-i(i+ft-i); 



Proof. Since the calculations are very similar to those in |BMRT10] . in order to 
avoid repetition, we will omit them here. □ 



2.2.4. Nature of critical points in various limits. Let us analyze the solutions to (18). 
Assume 5 > 0. The case 5 < is similar. Consider the two limits p — 00 and p — > in 
various scenarios depending on the angles and /3j. 



In the limit j9 — )■ 00 the solution to (|18|) has the asymptotics 
(19) s 
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When ]? — )• 0, 



(20) 



zbiTT 



p 



To get these asymptotics, it is necessary to show that if — < and p — )■ oo, 
or — > and p — )■ 0, then s ^ 1. This is a technicahty that has been addressed 
in jBMRTlOj . 



In Lemma [2.31 set 



c 

w 

Xi 

Vi 



1,2, 



1,2, 



/3 



..,n- j + 
= 1,2,...,J-1, 
1,2, ... ,n - J. 



Lemma 2.3 and its proof can be used to show that if z satisfies (15), then in the following 



situations it must be real. The idea is to show that if z is not real, then one of the angles 
in ^ with non-zero coefficient is much larger than all the other angles with non-zero 
coefficients, which is not possible. 



Case 1: p — )■ oo, s — )■ 0, f3j 
Also, 



and 



1. In Lemma 2.3 set ir = —1. Then (15) has the form of l\9h. 



IC 



w\ 



e'^lel + o{e) 



\w — Xi\ 



Hence, 



"'^ ^ — — s — )■ 0. The proof of Lemma 



e^\S\+o{S). 

implies that C must be real. 



2.3 



Since z = (, z must be real. 
Case 2: p — )■ 0, s — 7- oo, f3j = 1. The setup is similar to the first case and \( — w\ — »• 



\w 



Xi\ — > 0. Unlike the previous case. 



\w—xi I 



— oo. The proof of Lemma 



2.3 



les 



gives that in the limit considered /xi ^ a ^ (all the other angles), which imp 
( must be real. 

Case 3: p — )■ 0, s — )■ oo, jSj-i = —1, j3j ^ 1. In this case (15) has the form of ([9]) with the 
term roa in ^ replaced by — ^(1 + /3j)a + 1(1 — /3j)a. Notice that pi = 0. From 



the proof of Lemma |2.3| it is easy to see that if C is not real, then /^j — )■ for 

a ^ z/j for i > 1. Since 



i > 1, i/j — )■ for i > 1, a + a S> /ij for z > 1 and a - 
}'^~'^\ ~ ¥ = s — )■ oo, it follows that a — )■ (see Figure 



A calculation similar to ( 13 ) yields 



10). Thus, a — )■ as well. 



tan(a;) tan(Q;) fw-x\ + ^ 



tan(a) 



IC-^il 
\yi-x\\ 



+ 1 



s — )■ oo, 
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which imphes a ^ a. Thus, in the considered hmit a is much larger than all the 
other angles that appear with non-zero coefficients, which is impossible. Hence, ( 
must be real. 

Case 4: p — )■ oo, s — )■ 0,/3j = —1. Now (15) has the form of ([O]) if a is replaced by ce and 



tti = 1. In this case |C — w| — > 0, \xi — w| — )■ 0, and -^^z^ 0. It is easy to see 



from the proof of Lemma 2.3 that a ^ /ij, z/j and hence that C must be real. 



It is easy to see that the complex critical points Zcr of S{z) obtained from (19) and 



(20) are non-real complex except in the cases listed above. 



2.2.5. Critical points at finite (r, In the previous section the number of non-real com- 
plex critical points of Sr^^i^) identified when x is large. This section studies the 
number of non-real complex critical points for an arbitrary point (to,Xo)- Suppose that 
at this point the number of non-real complex critical points is neither 2 nor 0. It must 
be even, since they come in conjugate pairs. Assume there are 2m such critical points. 
The number of such critical points depends on the point (to,Xo)- However, if (ro,Xo) 
continuously changes in the r, x plane, the number of non-real complex critical points of 
Sro,xoi^) change, until it reaches a point were the equations 

(21) S'^Jz) = S':jz) = 

have a real solution 2; G M. In other words, the number of non-real complex critical points 
can change only near points (t, x) where Sr^xi^) has double real critical points. From ([s]) 
and ([7]) the condition (21) is equivalent to 

1 



(22) 
and 



Er=o^(A+i-A); 



X 



ze~^^-^ - 1 



i=i ^ ' 

(23) +^(l-/3.)lnn^Z73T +E^(l-/5^)l- 



2' ' " \ze-^-l ^2' ' \ze-^^-^-l 
-\r + \{V{V,) + V,). 

Think of this as a curve (r(z), x(^)) the r, x plane parametrized by 2; G M. Consider 
the complement of this curve in the r, x plane. The number of non-real complex critical 
points of Sr,x is the same for all points (r, x) inside the same connected component of 
this complement. Let D be the connected component which contains (tq, Xo) and let 
{t{zq),x{zq)) be a generic point on the boundary of D. The complement of (r(z),x(-2)) 
has another connected component whose boundary contains {t{zq)^ x{zo))- Call it D. The 
number of non-real complex critical points of 5*^-,;^; is 2m for points (r, x) £ 2) and 2m ± 2 
for points (t, x) G T). From Lemma 2.1 it follows that Zq must be in one of the intervals 
(e^'-i,e^^),/3i = ±1 or in (-00, e^^) U (e^",oo). Suppose, for example, zq G (e^'-i,e^'). 
The portion of the curve (r(z),x(^)) corresponding to z G (e^'-i,e^') is contained in 

2) n where D denotes the closure of D (a detailed study of the curve (r(z),x(^)) is 
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carried out in Section |4|). Now, when z approaches e^' (or e^'-^), x{z) will approach oo. 
However, it was shown in the previous section, that when x is very large, the number of 
complex critical points is either 2 or 0. Thus, 2m = or 2m = 2. 

This establishes the following proposition: 

Proposition 2.5. For any {t,x) the number of non-real complex critical points of Sr,x{z) 
is 2 or 0. Divide the t,x plane into regions according to the number of non-real complex 
critical points of Sr^x^z). Each connected component where the number of such critical 
points is has points where x is arbitrarily large. 



3. ASYMPTOTICS OF THE CORRELATION KERNEL 

This section analyzes the asymptotics of the correlation kernel Kx^^r.^{{t'l, h^), (tg, ^2)) 
in the scaling limit when Ymvk^oo^k = 0, the skew plane partitions are scaled by in all 
directions, 



4.k_+k „„J /\fh\ ■— hk Uk 



lim r^t^ = lim rkt2 = t, lim r^hi = lim r^/ig = X; 

A;— >oo A'— ^-oo A'— ^-oo fc— >oo 



and A(t) := ^1—^2 and A{h) := h\-~h\ are constants. A version of the saddle point method 
is used for calculating the asymptotics of the correlation kernel ([T]). The arguments used 
are along the lines of |UR03] . jpkoT] . |BMRT10] . 

Let r, X be such that S^-^x two non-real complex critical points. Deform the 
contours of integration Cz and Cw in the double integral representation of the correlation 
kernel given in @ioC'z,C'^in such a way that the new contours pass transversely through 
the two critical points of S{z) and ^{S{z)) < ^{S{w)) Wz G C'^^'im G C^, with equality if 
and only ii z = w = z^r- It was shown in the proof of Theorem 4.1 in |BMRT10] that this 
can be done for the function S{z) when /3j 7^ ±1,V2. The argument used there is general 
and applies for arbitrary /3j G [—1, 1]. During this contour deformation the contours cross 
each other along a path between the complex critical points of S'(z), so the residues from 
the term should be picked. Thus, the integral ([!]) can be written as 

^ i,k hk.. _ 1 f f '^-,b,,{z,t^^)<!>+,b,^{w,tl) 

K,,,,iit,,h,), (t„ h,)) - ^2^.^, y^^^^ J^^^^ X 



z — w zw 

fk\,f^ . +k\ 



+ 



^ r^^r, ^_^,^^{z,t'^)^^^,,^{z,tl] 
27riX (z,t^)<l>_,,, 



Recall that the first integral in the above formula has the form 

1 /" /" gT,x(z)-Sr,x("') I Q(^\ 1 

e '■fe dwdz. 



(2vri)2 Jc> Jc^ z-w 

S(z)-S(w) 

In the limit A; — ?■ 00 the asymptotically leading term of the double integral is e and 
3?(S'(z)) < ^{S{w)) along the contours C'^,C'^ except at the critical points. This implies 
that the main contribution to the integral comes from the critical points, but since the 
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contours of integration cross transversely at the critical points, the integral is zero in the 
limit. Hence, 



fe— ^oo 



1 /-^"^ (2; t^) ,fe ,fe 1, . fc. 1, , fe. 



Suppose tj' < Let Px^^it\,t2) = #{-D" n (^^,^2)}- The correlation kernel can be 
written as follows: 

fe— >oo 

= lim / — ; ; ^dz 

1 /"^ci^ 1 

= lim — / n (-e-^^^z) n , 



lim — / (-e-")^^fe(*^'*2)(i _e^"z)^*z-^'^+^^fe(*i'*')+5^'Afc(tJ)-|bAfc(t^)-i^^ 



fc— >-oo 27ri 



fc^oo 27ri 



In the case < ^i? K\k,qk exactly the same expression as above, with — -PAfc(^i,^2) 
instead of P\^(t\,t2)- It is easy to verify from the definition of -PAfc(^i,^2) ^^at it has 
the form f{t\) — f{t2)- Recall that the correlation functions are given by a determinant. 
Therefore, the process described by the kernel K\^^q^ is the same as the process with a 
kernel without the term (— e'^)^'^^'=*-*i'*2). Thus, the following theorem is obtained: 

Theorem 3.1. The correlation functions of the system near a point (r, x) 'the hulk are 
given by the incomplete beta kernel 

(24) K^A^t, Ah) = /(I - e-^^)^*z-^'^-^-i 



7 



27ri2; ' 



where the integration contour connects the two non-real complex critical points of Sr^z), 
passing through the real line in the interval (0, 1) if A{t) > and through (— oo,0) other- 
wise. 



When r, x are such that Sr,x has no non-real complex critical points, there are two 
possibilities. Either along the original contours Cz, it is true that Iki^Si^z)) < Ik^Siw)), 
in which case the above arguments imply the correlation kernel is zero, or otherwise, when 
the contours are deformed to obtain Iki^Si^z)) < 3fJ(5'(if)), one of the contours completely 
passes over the other and the correlation kernel is given again by (24) with 7 now being a 
full circle around and such that e"^ is not inside it. Taking Ah = At = gives that the 
limit of one point correlation functions, Ymvk^ao P\k,qki4^ ^1) either or 1. This means 
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that if Sr^xi^) non-real complex critical points, the tiles near the point t,x are 

horizontal tiles with probability or 1, and therefore the point (r, x) is in a frozen region. 



The above results imply that the curve {t{z),x{z)) where S'^ 



is 



the boundary between the frozen and liquid regions, so called frozen boundary. Section |4] 
studies this curve in detail. 



Theorem 1.2 follows immediately from Theorem 3.1, since (24) gives that the corre- 
lation kernel depends only on the critical points of the function Sr,x{^) which, as can be 
seen from its formula given in ([6]), is independent of the sequence of partitions and 
only depends on the function V{t) giving the back wall in the scaling limit. 



3.1. Correlation kernels on the frozen boundary and when x — )■ oo. The correla- 
tion kernel on the frozen boundary is given by the Airy kernel and near the cusps by the 
Pearcey process. The calculations are almost identical to those in |OR07j . therefore they 
will be omitted. 

It is proven in Section 4.2 of |BMRT10] that if in the limit x ~^ oo, S{z) has one 
pair of complex conjugate critical points and they have the asymptotics given in Lemma 
2.41 then the local statistics of the system is described by the bead process of Boutillier 



|Bou09j . 



The analysis in Section 2.2.4 shows that in the limit 



X oo, r = Vj-i + 6,6-^0, 

if 

p ^ gx-x^-i|^|i-5(ft-/3i-i) 

is fixed, then depending on the type of the angle at Vj_i and whether p ^ 1 or p ^ 1 
either the bead process or a frozen region is observed. This is illustrated in Figures [8] and 



4. The boundary of the limit shape 



The purpose of this section is to study the frozen boundary. It consists of those 
points (t, x) for which Sr^xi^) has double real critical points. In this section the most 
useful representation of the frozen boundary is as the parametric curve {t{z), x{^)) where 
t{z) and x{^) are given by (22) and (23), and 2; G M runs over all real numbers for which 
t{z) and x{^) are real. 



Define it to be 



il= ((-oo,e^°)U(e^",oo))|J |J (e^^S e^' 



VI <i<n 



It was established in Section 



2.2.2 



that if t{z) and xi^) are real, then z must be in it. 
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4.1. t{z) and xi^) ^ire real for all z E il. 



Proposition 4.1. For every z E ii, r = t(z) and x = xi^) given by (22) and (23) are 



real. In particular, Vz G it the point {t{z),x{z)) is on the frozen boundary. 

Before proving this proposition we need to prove one lemma. 
Define a function T{z) by 

n 

(25) nz):=J2o(^^-^^-^^)7^- 



2 

i=0 



t{z) from (22) can be expressed in terms of T{z) by 

r{z) 1 



(26) e- = .. 

Lemma 4.2. For any integer m > the 2m'th derivative of T , T^'^"^\z) ^ for all 
2 Gil. 



Proof. Formula (25) gives that 



T^'-\z) = J2li^^^l-P^) 



(_l)2m(2m)! 

i=0 ^ ^ 



Suppose by contradiction that 3zq G il such that T^'^^\zo) = 0. From the definition of 
il it follows that zq G {-oo,e^°) U (e^",cx)) or zq G (e^'-i,e^O with /3/ = ±1. The case 
zq G (e^'-i,e^'), /3z = 1 will be treated. Other cases are similar. 

When there are no corners, i.e. when n = 1, the proof is trivial. The rest of the proof 
is by induction on the number of corners of the limiting back wall. Suppose that there is 
only one corner, so n = 2, < < V2, Z = 1, e^" < 2; < (3i = 1, and /32 G [-1, 1). 
Then 



T(2™)(^o) = ((-l)'"^(2m)!) I 



/32 - 1 ^ 1 - /32 



2 \z - e'^o z- e^i z - e^^ 

^ ^2V-2-e^° (2-e^i)(2-e^2)y 

which contradicts the assumption T^'^'^^i^z^^ = 0. Thus, T^'^"^\zo) 7^ when the back wall 
has a single corner. 

Now, suppose the result is true when the number of corners is less than n — 1, and 
prove it when the number of corners is n — 1. First, consider the case n > / + 1, i.e. when 
there is at least one corner to the right of the corner at position VJ. It will be shown that 
if T(2™)(^q) = 0, the number of corners to the right of Zq can be reduced and still have 
that T(2™)(2o) = 0. 

Consider the function T{z) defined as 



2^' ^ ' 'z-e^- 2^ 'z-e'' 2'^ ^ ' z - e' 



6" 
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Notice that f{z,Vi,Vi+i) = T{z). In particular, g^T(^o, V^, = T(2™)(zo) = 0. 

Consider two cases. 

Case 1: /3i+2 — A+i > 0. In this case ■^^T{zq, a, b) is an increasing function of b for b 
such that > zq and thus -§^T{zo,Vi,Vi) < 0. From /3j = 1 it follows that 
|(/3i+2 — /?/)< which implies that ^^T{zo, a, a) is decreasing as a function of a 
when > zq and that 

lim ———T(zQ,a,a) = oo. 

Thus, 3ao G M such that Zq < < e^' and -^^T^Zq, oq, cto) = 0. This con- 
tradicts the inductive assumption since T( Z, flg, do) is equal to the function T{z) 
corresponding to a piecewise linear back wall with parameters Vq, Vi, . . . , V/_i, Oq, 
Vi+2, ■ ■ ■ yVn and /32, • • • , A; /9z+2, • • • , /3n, whlch has at most n — 2 corners (see 



Figure 12) 



Case 2: /3/4.2 — < 0. The argument is similar to the previous case. Now, -^^T{zq, a, b) 
is a decreasing function of b for b such that > and thus -^^T^Zq, Vi, V/_|_2) < 0. 
■^^T{zo, a, V1+2) is decreasing as a function of a when > 2:0 and 

a-s>ln{2o)+ aZ^™ 

Thus, 3ao G M such that zq < e"'^ < e^' and ■§^T{zq, oq, VJ4.2) = 0. This 
contradicts the inductive assumption since T{z, gq, V1+2) is equal to the function 
T{z) corresponding to a piecewise linear back wall with parameters Vq, Vi, . . . , V;_i, 
^o, Vi+2, ■ ■ ■ ,Vn and /3i, (32, ... , A+i? A+s? • • • ? /^n, which has at most n — 2 corners. 

This concludes the proof when there is at least one corner to the right of the corner 
at position V/. When there are corners to the left, the argument can be easily modified 
to show that if T^'^"^\zo) = 0, the number of corners to the left of zq can be reduced and 
still have that T(2™)(^p) =0. □ 
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Proof. [Proof of Proposition 4.1 Again, only the case z ^ {e 



with Pi 



4.2 



will be presented, as other cases work in th e sa me way. From (|26|), to show r is real it is 
enough to show that z — > 0. Lemma 

to T{z) 

nothing to show, so assume T{z) > 0. 



gives that T{z) ^ 0. If T(z) < 0, there is 



Again, do induction on the number of corners on the back wall. When there is only 



one corner, i.e. when n = 2, Vq < Vi < V2, I 
the inequality can be rewritten as follows: 



> ^ 1 < zT(z) ^l<-z 
T{z) ^ ' 2 



1, e^o < z < e^i, /3i = 1, and /^s e [-1, 1) 



+ 



/32 



- 1 1 

e^i z ■ 



/32 



3V2 



,V2 



< 



1. 



(2; - e^i)(z - e^i) z-e^^ 

The last inequality is true since the LHS is less than zero and the RHS is greater than 
zero. 



Now, assume that z 



7^ > whenever the back wall has at most n — 2 corners. Fix 

T{z) 

z and let < / < n — 1 be such that z < e'^K Such / exists if not all corners are to the 
left of z. The case when all corners are to the left of z can be treated similarly. 

Depending on the sign 



4.2 



Consider the function T(z, Vi, h) from the proof of Lemma 
of l3i+2 — A+I5 T{z, is either increasing or decreasing in 6 for 6 > VJ. He nce, either 
< f{z,VuVi) < T{z) or < f{z,Vi,Vi+2) 



< T(z). However, as in Lemma 



4.2 



both 



T{z, Vi, Vi) and T{z, Vi, are equal to the function T{z) corresponding to back walls 
with at most n — 2 corners. Thus, the induction hypothesis gives that z — 



1 



and z 



TR>0. 



> 



This establishes that if 2; G il, then t{z) = \n{z 



is real. 



Now, let us show that x{^) is real. This is obvious if z G 



z e ie 



V,. 



)U(e^", 00). Suppose 



with A = ±L If A = 1, (A = -1), then A+i - A < OjBi+i - A > 0). If 



e^'-e, it follows from (25) that T(^) > 0, {T{z) < 0). By Lemma 4.2, r(;z) 7^ 0, which 



implies T{z) > 0, {T{z) < 0) for all 2 G (e 



V,. 



Now, (26) implies that e'^ < z < e 



Vi 



(e' > z> e 



Vi- 



). Therefore, the coefficient of ln( 



in(23)isl-A = 0, (1 + A = 0). 



In conclusion, it was shown that if z G (e '~i,e ') with A = ±1, then the coefficient of 
ln( 



ze^l-l 



in (23) is always zero, which implies that is real for all z G it. 



□ 



4.2. The frozen boundary: connected components, cusps. 

Proposition 4.3. The frozen boundary consists of connected components, one for each 
outer corner where at least one of the slopes is ±1, and one connected component at the 
bottom (see Figure^. 



(27) 



Proof. If (r, x) is on the frozen boundary, then it is easy to see from (23) that 

dx ^ z _ 1 
dr ~ z - e'^ 2 
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Suppoze z G (e^'-i, e^') C it, for some /, for which f^i = ±1. Assume A = 1 (A = —1 can 
be done similarly). It is immediate from the formulas for t{z),x{z) that {t{z),x{z)) is a 
continuous curve when z G (e^'-^e^O- Consider two cases. 



First, suppose that ^ —1. Then (22) implies 
(28) lim r(^) = VJ-, and lim r(z) = Vi_i-, 



where for example by lim^_^gV;_ r(z) = VJ— we mean t{z) converges to VJ from below 
when z converges to e^' from below. 



Looking at (23) it is easy to see that 

lim x{z) = +00. 

z^e^i or e^i-i 

These calculations imply that when z ranges over {e^^-^,e^'), it gives a connected compo- 
nent of the curve {t{z),x{z)). 



If A-i = -1, then (l22j) gives that 

lim t(z) = lim t(z) = V/_i. 



Consider values of z which range over (e^'-^, e^'-^) U (e^'-^, e^'). It follows from (23) that 

/ _T/. 

Xo 



iim^xW = -i:^(i + ft)'"(;S^)- 

1=1 ^ ^ 

J\ 1 / - 1 \ 1 1 



i=/+l 

and 



lim x{z) = +00. 



z— ^-e^i or e^'-2 

These imply that when z ranges over (e^'-^, e^'-^) U {e^'-'^ , e^') , one connected com- 
ponent of the curve {t{z), xi^)) is obtained. The two pieces of this connected component 
corresponding to the intervals z G (e^'-^, e^'-^) and z G (e^'-i,e^') are connected at the 
point (V^_i,Xo)- 

Values z ^ [e^°,e^"] correspond to the component at the bottom. □ 



Figure 13. All possible corners. 




Next, it will be shown that each component corresponding to an outer corner develops 
a cusp. The above results, together with the analysis in Section 3.1 imply the frozen 
boundary at corners looks as described in Figure [13] 



24 SEVAK MKRTCHYAN 

A cusp appears if J = g = 0. This is equivalent to S'{z) = S"{z) = S"'{z) = 0. 
Lemma 4.4. Suppose (r, x) is on the frozen boundary. Then 

S"'{z) = ^ ^ = 0. 

dz 

Proof. By definition of T{z) 

1 



so 

dr T' ^, ^r, ^, 1 

0^1 + — = 0^T' + T2 = 0^0 = T' + 



dz {z-e^y 
However, notice that S"'{z) = T' + j^^. Thus, ^ = ^ S"'{z) =0. □ 

Theorem 4.5. Each connected component of the frozen boundary corresponding to an 
outer corner has a cusp on it (see Figures^ and 13). In particular, the number of cusps 
is equal to the number of outer corners where at least one of the slopes is a lattice slope. 

Proof. First, let us treat the case e^'-^ < z < e^', > —1, f3i = 1. 



By Lemma [4.4| cusps correspond to points z such that 

1 , dr 

and — = 0. 



T(z) dz 

This is equivalent to 

T(z) = and T'(z) 

z - \z -e'^ 

which means that to prove the theorem it is enough to show that when e^'-'^ < z < e^', 
then T{z) and -^—f are tangent to each other at a single point. 



Recall that T{z) is given by 

n 



2 

i=0 



It follows that \im^^^vi_i_^T{z) = oo = lim^_j^gVj_ T{z). Lemma 
0. In particular it follows that in the interval e^'-^ < 2; < e^' t 



4.2 



gives that T^^"'\z) ^ 
unction T{z) is never 



;ne 

zero, is always concave up, and has a positive minimum in that interval. 

When T > Vi, does not intersect T(z) in the interval e^'-^ < 2; < e^'. When r 
decreases, it intersects at a point. When r = V/_i, then 

lim T{z) = lim 77 — jr—^ + finite terms = 00, 

so and T{z) intersect. T{z) has a positive minimum in the interval e^'-^ < z < e^', 
so if T is small enough, -^—^ and T{z) will be tangent to each other. This completes the 
proof in this case. 
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Figure 14. 




The case e^'-^ < z < = — < 1 is identical. 

It remains to study the case e^'-^ < z < e^'-, (3i^i = — 1,A = 1. The same thing 
happens as above, except now the infinite terms in 

hm ^ T(z) 

cancel each other. If the limit is still positive, the argument works the same way as above, 
and a cusp r < V/_i is obtained. If the limit is negative, then as above it can be shown 
that a cusp r > Vi-i appears when e^'^^ < z < e^'~^ If the limit is zero, the cusp is at 
r = Vi-i. □ 



Notice that (28) implies ^ — ?■ ±00 when z — )■ e^',VZ. Examining (27) it can also be 



estabhshed that ^ 7^ 0, and since for z G (e^'-i,e^'), ^^^^ is a continuous function, it 
follows that in the intervals z G (e^''i,e^'), ^^^^ does not change its sign. 

In the r, x plane each connected component except the one at the bottom starts at 
T = Vi-i,x = 00 for some i, x decreases until it reaches the cusp, then increases to 



Vi,x = 00. See Figures P and 13 




The results of this and previous sections give the following characterization of the 
frozen boundary: 
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(1) On the pieces of the back wall where the slope is ±1 the disordered region is 
bounded above, while at other places it grows infinitely high. 

(2) The number of connected components of the frozen boundary is one more than 
the number of outer corners where at least one of the slopes at the corner is a 
lattice slope. 

(3) The frozen boundary develops a cusp for each such outer corner. 



Figures [9, 15 and 16 give several interesting examples of frozen boundaries. 



Figure 16. V = 0.7, 1, 1.05, 1.2, 1.25, 1.5, 1.55, 3, 3.1, /3 = 1, 0.7, 1, 0.7, 1, 0.7, 1, -1. 




Appendix A. An integral formula for ^^^(-z) 

One of the main difficulties in this paper, and in fact also in jOR03] and |OR07j . is 
understanding the asymptotics of ([T]). Similarly as in |OR03] we study the asymptotics 
of the function 

Tk In $b (z, tk) = Tfc In ' ' 



in the limit when — )■ and tk is a sequence such that rktk t (here ^-,6^^. i^^ ^) ^^"^ 
t) are as in (|2|). In this appendix we prove a technical result concerning this 
asymptotics in a very general setting. We consider an arbitrary continuous function 

V{t) : M ^ M, 

which is Lipshitz with constant 1. This is a less restrictive assumption on V than in 
the rest of the paper, where V{t) is assumed to be piecewise linear, and the domain 
is restricted to a finite interval. Otherwise, we use notation from the introduction. In 
particular, we are considering a sequence of back walls bx^, {t) and G ]R>o such that 
tends to 0, and functions Bx^{t) defined by 



converge uniformly to V{t). 
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Lemma A.l. If Bx^. uniformly converge to V , liiiifc^oo ''^fc = 0, andlivuk^oofktk = then 

_ -,M^-1 

r,ln$ft, (zAi.) = UT + V{T))\n{l-e^z-^)- / Um + V{M)) — 



1 r I -e^^z~^ 

lim r. In {z, t,) = -(r + \/(r)) In (l - e^z-') - j -(M + V{M)),^—^^dM 

+ -(r-\/(r))ln(l-e--z)+ y -(M - T/(M))^— ^^dM. 

/n j)artzc«/ar limfc_!.oo In $5^^ (z, tfc) is independent of the family {h\^{t)}k. The conver- 
gence is uniform in z on compact subsets of C 

Proof. Let us analyze the denominator of t^). The same arguments will 

work for the numerator. Define 

Dk := rk In 

From Q we get 



^+,bAz,tk) 



Du = -rfc I In I n (1 - = - E (1 - ^"'■*"^) 1 • 

ym>tk,'meD+,meZ+l^ / / \ m>tfc ,6'^ (m)=-l 



Notice that 



Make the change of variable M = rkm and set = rk{tk + \)- Dk can be rewritten as 

Dk = -Yl ^4(^ - ^'^•("')) - ""'"^) = 



'2 



5^ rS{l-B',{M))\n{l-e-''z). 



'2 

Me {rt: ,Tfe +rfc ,Tfe +2rfc , . . . } 



To make formulas simpler, define 

/(M) := - In (1 - e-^^^) . 

Let a > r, and let cxfc be a sequence such that ^^7^ £ ^+ and limfc_j.oo cta: = cr. Define 

^.^:= E r,i(l-i?UM))/(M). 

Since Bxi_{s) has constant slope when s G (M, M + r^) we can write D'^ as 

= ^ E ([M + r, - S,,(M + r,)] - [M - B,,(M)])/(M). 

Me{Tfc,rfc+rfc,Tfc+2rfc,...,o-fe} 
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Now, using {ttn+i - an)b„ = (a„+i6„+i - a„fo„) - a„+i(&„+i - &„) rewrite D'^ as 
Dl^l m+n-B,,(M+ n))f(M + n) - (M - B,, (M))f(M)] 

■We {rfc ,Tfc +rk ,Tk+2rk,-,<7k} 

-I E (M + rfc - S,,(M + r,))(/(M + r,) - /(M)). 

Me{Tk,Tk +rk ,Tk+2rk,-,Tk} 

Summing up the telescoping sum and rewriting the second one give 

= ^{(^k + rk- Bx^{(Tk + rk))f{(Tk + Tk) - ^(Tfc - Bx^{Tk))f{Tk) 



1 

2 



rk{M + rk- Bx^{M + rk)) \^ — j. 



M&{Tk,Tk+rk,Tk+2rk,---,crk} 

Since B\^{M) converge uniformly to V{M), and f'{M) and f"{M) are bounded when 
M e (rfe, (Tfe), we have 



lim 



^fc - ( li^k - V{ak))f{ak) - liTk - V{Tk))f{rk) 



\ rk{M-V{M))nM) 

Me{Tk,Tk+rk,Tk+2rk,-;(Tk} 



0. 



The sum above is a Riemann sum for | (M — V{M))f'{M)dM and thus in the limit 
rjfc ^ it converges to the integral, so we have 



Since limg-^oo ~ = 0, we obtain 

^lim = (-^(r - V{T))f{T) - £ \{M - V{M))f{M)dM^ . 



Similarly, if A^^. is defined as 

Nk := Vk (in , 

it can be shown that 

^lim AT, = +-(r + V{t)) In (l - e^^-^) ' j ^^(M + V{M))^—^^dM. 



Combining the results completes the proof. □ 
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Corollary A. 2. IfVij) is piecewise differentiable, liiiifc^oo ''^fc = 0, and limk^oo f^ktk = t, 
then 

lim rfcln$fe, = - [ -^(1 + V (M)) In (l - e^V^) 

POO 1 

+ / --{l-V'{M))ln{l-e-^'z)dM. 

Jt 2 



Proof. This follows from Lemma A.l by integration by parts. □ 
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